Calculating the Fine Structure of a 
Fabry-Perot Resonator using Spheroidal 
Wave Functions 

M. Zeppenfeld' and P.W.H. Pinkse^ ^ 

^ Max-Planck-Institut fiir Quantenoptik, Hans-Kopfermann-Str. 1, D-85748 Garching, 

Germany 

^Mesa+ Institute for Nanotechnology, University of Twente, PO Box 217, 7500 AE Enschede, 

The Netherlands 
email: martin. zeppenfeld@mpq. mpg. de 



Abstract: A new set of vector solutions to Maxwell's equations based 
on solutions to the wave equation in spheroidal coordinates allows laser 
beams to be described beyond the paraxial approximation. Using these 
solutions allows us to calculate the complete first-order corrections in the 
short-wavelength limit to eigenmodes and eigenfrequencies in a Fabry- 
Perot resonator with perfectly conducting mirrors. Experimentally relevant 
effects are predicted. Modes which are degenerate according to the paraxial 
approximation are split according to their total angular momentum. This 
includes a splitting due to coupling between orbital angular momentum and 
spin angular momentum. 
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1. Introduction 



From cavity quantum electrodynamics to deterministic single-photon generation, a host of ex- 
citing recent experiments are made possible due to the precise control of light achievable inside 
optical resonators {HO. In the case of Fabry-Perot resonators, the detailed understanding of the 
resonator modes necessary for these experiments is generally based on solutions to the paraxial 
wave equation. Despite being approximate in nature, these solutions allow many crucial fea- 
tures of the resonator modes to be explained, from the Gaussian profile of the lowest-order 
TEMo.o mode to the regularly spaced eigenmode spectrum. 

Despite its success, novel experiments allow the boundary of the validity of the paraxial 
approximation to be probed. This includes microresonators, where the small mode volume 
leads to a violation of the assumed paraxiality, and ultra-high-finesse resonators, where the 
high spectral resolution allows even minute details to be resolved. Understanding all observable 
features in such experiments requires the electromagnetic field to be treated more precisely than 
the paraxial approximation allows. 

In the past, several approaches have been used to obtain corrections to the paraxial approxi- 
mation for resonators |[3l|4l|5]|6l. Lazutkin ||3| chose an analytic expansion of the wavefunction 
to satisfy the wave equation with appropriate boundary conditions. Laabs et al. H calculated 
a resonator round-trip propagation matrix using the method of Lax et al. Ill of calculating cor- 
rections to the paraxial equation. The effect of spherical aberrations was calculated by Visser et 
al. m based on perturbation to Gaussian propagation inside a resonator Zomer et al. ^ used 
a diffraction integral to propagate light around a resonator 

None of these papers take into account all effects needed to calculate the complete first-order 
correction beyond the paraxial approximation of the electromagnetic eigenfrequencies of a res- 
onator Lazutkin and Zomer et al. restricted their analysis to scalar fields in two-dimensional 
resonators. Visser et al. focused solely on spherical aberrations, albeit leading to an analytical 
result reproduced here. The most rigorous approach so far by Laabs et al. used an insufficient 
approximation by treating the mirror surfaces purely in terms of a position-dependent phase 
shift. 

An alternative approach for treatment of the resonator eigenmode problem is through the use 
of spheroidal coordinates. The spheroidal coordinate system, depicted in Fig.[Tl is ideally suited 
for calculations with Gaussian beams due to the match between surfaces of constant phase of 
a beam and surfaces of constant £, as well as between the variation in beam diameter along 
a beam and surfaces of constant 77 . With appropriate approximations, these coordinates have 
been applied in the past to resonators to obtain results in agreement with paraxial theory miH. 

In this paper, spheroidal coordinates are used to calculate the eigenfrequencies of a Fabry- 
Perot resonator to first order beyond the paraxial approximation in the short-wavelength limit. 
The electric field inside the resonator is expanded in terms of vector spheroidal wave functions. 
Requiring the boundary conditions for perfectly conducting mirrors to be satisfied allows the 
expansion coefficients of the outgoing wave at each mirror to be expressed in terms of the ex- 
pansion coefficients of the incoming wave. The resulting round-trip propagation matrix for the 
expansion coefficients is solved via perturbation theory, resulting in a compact expression for 
the resonator round-trip phase shift, Eqs. ( l40l i and (|4TI) . The paper concludes with a discussion 
on the experimental implications of this result. 

2. Mathematical Foundations 

This paper relies heavily on the theory of spheroidal wave functions as found elsewhere ifTOl 
fm . In particular, the theory of spheroidal wave functions in the short-wavelength limit has re- 
cently been significantly expanded in Ref. ifTTI . Here we briefly summarize the most important 
definitions and results from this paper. 




Ti=o: 



n=o.3 



Fig. 1. The oblate spheroidal coordinate system. The (^-coordinate describes the set of 
ellipses with a common pair of focal points separated by a distance d, the Tj -coordinate 
describes the set of hyperbolas with the same two focal points, and the (^-coordinate de- 
scribes the rotational angle around the z-axis. A pair of mirrors forming a resonators can be 
matched to the coordinate system as indicated, in this case with = — 1 .25 and = 1 .75 



The mapping from spheroidal coordinates , T] , ) to cylindrical coordinates (r, z, ) is given 

by 

r = \^J{^-n^){^ + ^^) and z = (1) 

the -coordinate being the same in both coordinate systems and d being the interfocal distance 
as shown in Fig. [T] The scalar wave equation is separable in spheroidal coordinates, allowing 
the solutions to be written as the product of three functions depending only on ^, Tj and ^, 
respectively. This leads to so-called scalar spheroidal wave functions ^,„y = /?,„v )Sn,v {T])e""'^ 
with 

and 

labeled by the indices m and v. Note that m is simply the integer orbital angular momentum due 
to the e'""'' ^-dependence. The variable x = 2c(l — f]), not to be confused with the Cartesian 
coordinate, is introduced to simplify calculations. The parameter c = kd/2, with wavevector k, 
quantifies the scaling of the coordinate system relative to the wavelength. For short wavelengths 
compared to d, c is a large number, and the functions r,„v ) and i,„v [x) can be expanded as 
asymptotic series in |. To lowest order in 4, r„,v(^) — 1 and Si„y{x) ~ Ly"^ (x) where Ly"' (x) is 
a Laguerre polynomial. 

Based on the scalar spheroidal wave functions, vector spheroidal wave functions E^^^, prop- 
agating in the ±£, direction are defined which satisfy the wave equation as well as V • Ey^^, = 0. 
The Ej^y can therefore be regarded as the electric field of a solution to Maxwell's equations in 
free space. With spin angular momentum a ~ a* corresponding to left and right circular polar- 
ization, y = m ± 1 denotes the total integer angular momentum of the field about the symmetry 
axis. 



3. Satisfying the Boundary Conditions 

The electric field inside the resonator is expanded in terms of vector spheroidal wave functions 
as 

E = L (biav^iav + bJav^Jav), (4) 
Jav 

with expansion coefficients bf^^. For perfectly conducting mirrors, the component of the elec- 
tric field parallel to a mirror surface S must vanish at S, i.e. Ejg y = 0. 

The mirror surfaces are matched to the spheroidal coordinate system as exemplified in Fig.[T] 
as follows. For a resonator of length L with spherical mirrors with radius of curvature R 
and R+, the spheroidal coordinate system is appropriately scaled by choosing the distance d 
between the focal points to be |8] 



4L{R++R^-L){R+-L){R^-L) 

"^-V {rTTW—W • (5) 

In this case, the two surfaces of constant £, = ^± with 

^^^^d{R++R^-2L) 

are separated by a distance L and have a radius of curvature of R± on the resonator axis. 

For a spherical mirror, a correction to the spheroidal surface of constant ^ is generally nec- 
essary. Switching briefly to cylindrical coordinates (z,r, ^), the mirror surfaces are specified as 
the distance z ~z±{r,^) above the plane z = 0. For a smooth, cylindrically symmetric mirror, 
z± ('", ^) can be expanded in powers of as 

z±(r,(^))=z±Tz^TC4±r4 + ^('-')- (7) 

ZK± 

Here, z± and C4± are expansion coefficients and R± is the radius of curvature of the surface at 
r = 0. For a spherical surface, C4± = 1/ (8/?|). Other values of C4± can be chosen to describe 
cylindrically-symmetric mirror aberrations. For example, for a spheroidal surface of constant 
^ = ^± we have z± = d^±/2 and C4± = l/{4d^±R\). For a parabolic mirror, C4± = by 
definition. 

Switching back to spheroidal coordinates, the mirror surfaces can equivalently be specified 
as the value of ^ as a function of x = 2c(l — Tj) and (j). Eq. ^ is transformed into 

2 / 3 



with 



^^ = ^±(x,.^))=^±T/4±^ + ^( ^ ) (8) 

d^lRl t± 
U±^^Y^C4±Ty^ (9) 

S± being the mirror surface at Note that for small r, x is proportional to cr^. Additionally, 
due to the factor in the definition of S„,v{t]), Eq. (O, the spheroidal wave functions vanish 
for X 3> 1, motivating the use of x in Eq. ([8]l. The term linear in x is missing due to the specific 
choice of d and ^± in Eqs. ^ and (|6]l. As will be seen, the x^ term in Eq. ^ is the highest- 
order term in x which must be retained to calculate first-order corrections to the resonator 
eigenfrequencies. 



The mirror surfaces can now be parameterized in terms of the transverse coordinates Tj and 
^ as 

(10) 



S±(tj,0) =x(|±(x(tj),0),tj,0 



where x((^ , 77 , 0) is the vector in M? from the origin to the point denoted in spheroidal coordi- 
nates by , T], 0). A basis for the tangent space to the mirror surface at a point on the mirror 
surface, needed to calculate the component of E parallel to the mirror surface, is given by 



(11) 



where the e„. are unit vectors and the /?„. are scale factors in spheroidal coordinates |[T2l . The 
component of the electric field parallel to the mirror surface is proportional to the inner product 
of the electric field with these basis vectors. 

We now show that as long as we are only interested in corrections to the resonator eigen- 
frequencies to first order in i, we can use as an approximation for the basis (fTTI) . For 

x= 1^(1), we have the following order of magnitudes. 



ax 



dx 



hi, 



(12) 



dn] 



as can be seen from the corresponding definitions in Ref. ifTTl . As a result, the contribution to 
• E of the • E term is of order c^^ compared to the contribution of the §7] • E term and can 
therefore be neglected, justifying our simplified basis. 

To calculate Ct] • E^^^^, and • E^^^^, we additionally need the following expressions. 



• (x±;y) = ±ie 
e^-(x±/y) = -77 

e^-(x±iy) = ^ J 



±1(4 



• Z 



/ 1+^2 



1 - „±<* 



e^-z = Tj^ 





772 


772- 








772- 





(13) 



For first-order corrections, the expressions with are not needed, but have been included for 
completeness. Using previous expressions and with a significant amount of algebra, one obtains 
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7(7+ V 



S± 



s± 



1-^ 

8c 



-r f ^2r W I f 2riJ) 



/(v + y-l)(2v-l) ,(j) ^ /(2v2 + 2/v-5v-2/+l) ^(j) 



(14) 



v-r 



i(2v2 + 27v + 5v + 37+ /(v+ l)(2v + 2y+ 1) ,(y) 



"v+l 







S± 


j- 




Jx 




1- 


8^ 



s± 



AJ) 



.(J) 



T/4±x2(V + /+ 1)4-" ±/4±x2(v+ I)l(^) 



/(v + y)(v + /+l)(2v + l) ,7) 



v+1 



(15) 



/(v + /+l)(2v^ + 3v + / + 2) ^(y) /(v + l)(2v^ + 47v + 5v + 27^ + 67 + 4) ^(j) 
8 + c ^v+r 



/(v+l)(v + 2)(2v + 2/ + 3) ,(7) ■ 

o ^v+2 



s± 



-E 



7CT+V 



S± 



Jx 



-x/2ij^_ 



,(J) 



1+a 



^y.v-i(^±)((v+/-l)L! 



iJ) 



,(/) 
+ 1 



/^y- 1 , V ) ( V + / - 1 )l[,-'12 - (3 V + 2/ - 1 )Lt-'Z 1 + (3 V + / + 1 )Li;' - ( V + 1 )L' 



(^) 



.(7) 



2(1+^1) 



and 



■(]) 
v+l 



(16) 



Ja^v 



8c,' 



/?yv(^±)(v + /+l) -(v + /)4-'Zi + (2v + /+1)4-''-(v+1)l: 



Rj,y+,{^±){v+\) ( - (v+y+ 1)4-" + (2v + y + 3)L<-'|i - (v + 2)L«2 



1+^1 

(v+l)(3v + 2y + 4)4-'ji + (v+l)(v + 2)4^2 



"v+l 



j-Rj+xy{^±) -(v + y)(v+/+l)Lfli + (v + /+l)(3v + / + 2)Li;'- 



■ iJ) 



2(1+^1) 



(17) 



The Ly''' are Laguerre polynomials L^/' [x) with the dependence on x being implicit. Heavy 
use of recursion relations among the Laguerre polynomials has been made. Note that we have 
evaluated the spheroidal wave functions at the mirror surfaces S±. For Ej^^y, the field running 
in the other direction, the results are the same except that ^± is replaced by — ^± and f^± is 
replaced by — /4±. 

Combining Eqs. (I14m7b with Eq. (HI, the boundary condition for E for the mirror at ^± can 
be reformulated as a matrix equation for the coefficients bj^^ of the outgoing wave in terms of 
the coefficients b^^^ of the incoming wave. We begin by defining the set of vector functions 
\jav as 



Jx 

si 



iJ) 

v-l 



(18) 



and 



yjo-v=^'^ (l - I) e-^^' ((V +/+ 1)4-" - (v+ (-/e^ - e,; 

and by defining the functions Uiy^v and U2jav according to 

\Rj.v-i{^) yjG- ,v-2 - "^Ja-.y- 1 - ( V + / - 1 )yja+.v- 1 + VV^^- 



1 



Rjv ) ( Vya- ,v- 1 - Vya- v - ( V + ^) Vya+ v + ( V + 1 ) Vy 



'o+,v+l I J 



(19) 



(20) 



Uiy<T-v = \/^((v + /+l)v/iv + (v + l)v/j,v+ije,, = 

i 7?yv ) ( V + 7 + 1 ) ^Vj„- , V- 1 - Vya- V - ( V + Vy^+ v + ( V + 1 ) Vyc,+ , v+ 1 
^ ^y,v+l ) ( V + 1 ) ( Vy^- V - yja-y+l - ( V + / + 1 ) Vy„+,v+ 1 + ( V + 2)Vy^4 



(21) 



U2ya+v =-'LTy-i,ve"''e,, = 

X 1 /1\ (22) 

- ^^y-l,v(^)Vya+v - ^^y-l,v(^)(Vya-,v-2-2Vy^-,v-l +Vy^-v) + ^ - ) 



and 



U2y<7-v =^V'i+i,ve" ''e,, = 
1 



^TJy+i ,v )Vy„- V - ^^7+1 ,v ) [( V + /) ( V + / + 1 )Vy„ 



(23) 



2(V+ 1)(V + /+ l)Vy„+ v+l + (V+ l)(v + 2)vy„+ ,+2] + ^ (i) ' 



The parallel component of the electric field of Ej^.^ at the mirror surfaces can then be written 
in terms of the v y^v as 



E 



7(7+ V 



It 



8c 



Vya+v +^y-i,v^/c7+,v-l +^y-i,v^/cr+,v+l 



1+^1 



= Uiy(j+v 



(24) 



and 



EycT-v 



= ^y+i,v(^±) 



It 



1+^ 



=Uiya-v 



1 ^1 



c 8c 
U2y 



Vyfj-v +^y+i,v^y(j-,v-i +^y+i,v^y(j",v+i 



^=^± 2c-i+^r2^^"^ 



(25) 



For E,j^^, we must again replace ^± everywhere by — ^± and f4± by — /4±. The Aj^ are expan 



sion coefficients for the angular spheroidal functions ifTTI . We have Aj^j ^ = ^i];) so that to 
lowest order in 4, Ey^^^, jg^ y is simply proportional to \jav, motivating the choice of \jav 
Eqs. ( l24l i and dZSl l have the form 



E 



7(7 V 



±l,±2J,a',CT 



(26) 



a'v' 



with appropriate coefficients The subscripts on the ± signs denote two indepen- 



dent choices of + or — . Eqs. ( l24l i and ( |25] l contain factors of x as part of the coefficients of the 
functions Vjav- These can be removed using the relations 



-^•V7CT+V = -(V+/-1)V, 



+ (2V +/)Vj„+v - (V+ l)Vy^+ v+1 



and 



{v+J+l)yj^ 



.v-l 



- (2v+/ + 2)vy^-^, - (v+ l)vy^-,v+i, 



(27) 
(28) 



making the coefficients a^,^^^^'''" independent of x and (j). 
The boundary condition for E can now be written as 



E 



s±,|| 



S±, 



' ^Jav^Jav 



a'v' 



L( +.±J.a'.a,+ 



S±, 



-.±J,a'.a,- 



(29) 



0. 



Due to cylindrical symmetry, modes of different J do not couple, so we have restricted our 
attention to a single /. The Vy^v, considered as vector functions of x and 0, are linearly in- 
dependent. As a result, the expression in the square bracket in Eq. ( |29] ) must be zero for each 
value of a' and v'. 

We define the following matrices and vectors. 



^a'a ~ y"v'v )v'v, 



1±1.±2 _ 



l±l.±2 



1±1,±2 



and 



(30) 

The matrix A*i'*2 maps the traveling wave described by the vector b(j' of coefficients of wave- 
functions Ej^,, onto a set of coefficients of surface functions v jav describing the electric field 
component parallel to the surface at S±2 . This finally allows us to write the boundary conditions 
as matrix equations. 



i-,± 



= 0. 



(31) 



Solving for bj in terms using the boundary condition imposed by the mirror at we find 
the round-trip matrix for the resonator to be given by 



(A' 



(32) 



The eigenvectors of A are vectors of coefficients bj^^ of resonator eigenmodes, the cor- 
responding eigenvalues are the round-trip phase shifts. To lowest order in 4, the matrix A is 
diagonal. The resonator eigenmodes are therefore of the form ^jfj±v^/(j±v ^ ^ja^v^ja^v '^^^ 
the corresponding round-trip phase shift is equal to 



^ 2/(f(,5+-^_)-(2v+(yTl) + l)(arctan(§+)-arctan(,?_)) + ^?'(i)) 



(33) 



This result is equivalent to the one obtained by considering a resonator in the framework of the 
paraxial approximation. 



4. First-Order Corrections to tlie Round-Trip Pliase Sliift 

First-order corrections to the round-trip phase shift of Eq. (l33T l can be obtained via degenerate 
perturbation theory. We begin by determining those eigenmodes of the round-trip matrix A 
which are degenerate to lowest order in 4. As can be seen from Eq. ( |33] |, the lowest-order 
eigenvalues corresponding to the modes Ej^j+y Ey^j- are the same for arbitrary J and 
V. This reflects the fact that, within the paraxial approximation, the resonator eigenfrequency is 
independent of the polarization. Additionally, for a cavity geometry such that ^"^'™('^i)^~^'''='™(^2) 
is a rational number, say ^ with relatively prime integers p and «, the lowest-order eigenvalues 
corresponding to the modes Ejav and Ey^ v+n are the same for arbitrary J, v and a. 

Lowest-order degeneracy between modes is only relevant if the modes are coupled by first- 
order matrix elements. As can be seen from Eqs. (l24l i and dZSl l together with Eqs. (l20H23b . a 
mode Ejfj+y in general couples to the modes Eyfj+ yj_i as well as to the modes Ej^j-yi with 
v' = V — 2, V — 1 and v. A mode E j„- y couples to the modes Ej^- ,v±i as wefl as to the modes 
Ey<T+v' with v' = V, V + 1 and V + 2. For nonzero f4± coefficient, the modes Ey^^+y and Ey^j-y 
additionally couple to the modes Eyo.+ y_|_2 and Ej^-y^2^ respectively. 

From this analysis, we see that the degeneracy between the modes Ey^^+y and Ey^j- y_j 
must be taken into account. On the other hand, a degeneracy between the modes Ey^^+y and 
E/(T+ v+n '^^^ to resonator geometry is only relevant for « = 1 and for « = 2, f4± ^ 0. For these 
two special cases, an infinite set of degenerate modes are coupled by first-order off-diagonal 
matrix elements, and finding the first-order corrections to the eigenvalues of A is significantly 
more difficult. We therefore exclude these two cases from further analysis. 

Disregarding the two cases « = 1 and n = 2, f4± ^ 0, the first-order corrections to the eigen- 
values of A can be obtained from the two-by-two submatrix of A which corresponds to the 
subspace spanned by E jfy+ y and Ey^j- y _ i . We denote this matrix as the two-by-two round-trip 
matrix. Since the four matrices A^i'*2 are all diagonal to lowest order in 4, the two-by-two 
round-trip matrix is given to first order in 4 by a product as in Eq. ( |32| |. except that each of the 
matrices A^' '^^ replaced by an appropriate two-by-two submatrix. Specifically, these four 
submatrices are made up of the coefficients a^^' '*^ with i,j G {0, 1}, given by 

±1.±2 _ ±1. ±7,7,(7+, (7+ ±1,±2 _ ± I, ±2 ,7,(7^,(7+ 

ago — flvv J «io —^v-l.v ' ^2^^ 

±l,±? ±i.±2,y,(7+,(7^ ±i,±2 ±|,±2,7,a^,a^ 

^01 —"v.v-l ' "^11 —"v-l.v-l 

The coefficients ajj''^, obtained by applying Eqs. dZTI ) and ( |28] ) to Eqs. (l24l) and ( 1251 ) and com- 
paring with Eq. ( l26l ). are given explicitly by 

+,± „ ,^,_ //4±(6v' + 67v + /^ + /) , 2v+J 



^/?y,y_l(^±)-^7;yy(^±) 
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01 
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4c 1+^2 



+ l.v-l(^± 



v(v+7) ^ 



2c 1+^ 



-^y+i.v-iCi^i) - 



For a^j"^ one can check that a^y* = [ajj^)* with * denoting complex conjugation. 

The half of the two-by-two round-trip matrix corresponding to the mirror at £,± is given by 



(35) 



no 



^01 



ho 
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=F.± T.± 



=F.± ±.±_ 

±.± 
'^10 '^00 



^00 '^10 



=F.± ±.± 
'^11 '^0^_j_' 



2jQ AqJ 



(36) 



Surprisingly, the off-diagonal elements in Eq. ( [36] l are zero to first order in \. This is due to 
the fact that the complex phase of the three functions i{R j y- 1 ) ^ ^/v )), ^7+i,v- 1 ) and 
Rj-i,v{^) is equal to lowest order in 4, i.e. 

arg(i(/J,,v-i(^)-/?jv(^)))=arg(/;y+i,v-i(^)) + ^('^) = arg(«y_i,v(^)) + ^ f ^ V (37) 



This in turn results in equal complex phase to lowest order in i for the four coefficients a^^ 
with fixed ±i and ±2. Since a^'* is the complex conjugate of fl^'*, the lowest-order term of 



±i,±i 



,± +,± 
fl., ., 
'J 

i 7^ we have 



a-' a-,'., is real for arbitrary j, i' and /. Since a-j' = for / — j and a-j'^ = ^(4) for 



00 "^lO 



'10 "00 



T,± ±,± 
'^00 '^10 



(fl 



■ a, 



00 "^10 



=F.± ±.± 

2qQ fljQ 



(38) 



and similarly 



fljj flQJ 



-flgj fljj 



For the diagonal elements in Eq. ( 1361 ), we note that flQj''*^fl]'^Q = so the top left ele- 

ment, which we denote is given by a^'^ /a^'^, and the bottom right element, which we 

denote is given by flj''j'''/fl^j'''. 

The fact that the off-diagonal elements in Eq. ( |36] | are zero through first order in 4 means 
that the mirrors do not couple the modes ^ja+v and Ey^^- y_j to first order. As a result, Ey^y+y 
and Ey(j- y_ 1 remain lowest-order resonator eigenmodes when taking into account first-order 
corrections, despite being degenerate at lowest order The round-trip phase shifts for the two 
modes are given by the products Aj„+y =A+^+^A7^+^ wAXj„-,^_i=X+^_^^Xj^_^_y With 



(39) 



a bit of algebra we obtain the central result, 



exp 
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(40) 



and 



exp 



±2iyc^±-{2v+J)civaeLn{^±)- 
(v+y)±/4±(6v(v + /)+/(/-l))^+^(^^^^ 



(41) 



Note that the first line of both Eqs. (|4Q] i and i4T[ is simply the paraxial result of Eq. ( [33] l. 
The term proportional to f4± has been calculated previously 13], and is simply equal to the 
phase shift corresponding to the intensity weighted average position shift of the mirror due to a 
change in f4± . 

5. Discussion of the Results 

The resonator round-trip phase shift of Eqs. ( l40l i and i4T[ allows the calculation of the exact 
next-order correction to the eigenfrequencies of a Fabry-Perot resonator with perfectly con- 
ducting mirrors. Beyond providing a bound for the accuracy of the paraxial approximation, 
the correction term has a number of physical implications. One of the most interesting is the 
prediction of spin-orbit coupling in a resonator For two modes with equal v and equal orbital 
angular momentum m = J^l but opposite spin a*, we find a difference in round-trip phase 
shift equal to 



nj+l,fJ+V 



exp 



±- 



exp 



2im \ 
kR^J 



(42) 



The resonator eigenfrequency is therefore equal to the sum of a term independent of spin and 
a term proportional to the product of spin and orbital angular momentum, i.e. spin and orbital 
angular momentum are coupled. 

The correction term for the resonator round-trip phase shift fundamentally alters the no- 
tion of degeneracy of modes in a resonator. In the absence of this term, all modes with equal 
2v + m are degenerate, independent of polarization. This allows significant arbitrariness in the 
choice of eigenmode basis for a resonator, i.e. Hermite-Gaussian modes, Laguerre-Gaussian 
modes and many other linear combinations are equally valid. For a resonator geometry such 
that £^ additional degeneracy results, making, e.g., a confocal resonator an 

apparently ideal choice in applications where a cavity with many degenerate modes is desired. 

Taking the corrections into account, this picture breaks down completely. Degeneracy be- 
tween the various modes is almost completely lifted. Circularly polarized Laguerre-Gaussian 
modes are clearly singled out as the fundamental resonator modes, consistent with the cylin- 
drical symmetry. This removal of degeneracy is the ultimate situation for spherical mirror 
resonators. However, the presence of the /4± term in the round-trip phase shift opens com- 
pletely new possibilities via the capability to manufacture aspherical mirrors. For a choice of 
/4± = =F^(^±/(1 the correction term for the round-trip phase shift is significantly re- 

duced, leaving only a term which is independent of V and mainly quadratic in J. This allows 
for the construction of resonators with a degeneracy of a large number of modes reestablished. 
Note however that this possibility does not include confocal resonators. Such resonators fall 



in the class of resonators with n = 1 as defined above, causing the derivation of Eqs. (l40l i and 
dH) in section m to fail. Together with the fact that they lie at the edge of the zone of stability, 
confocal resonators turn out to be among the worst possible choices for a degenerate resonator 
Last but not least, we consider the prospect of observing the effects of the correction term 
experimentally. While this should be relatively easy for a resonator with dimensions on the or- 
der of the wavelength, an alternative highly interesting system would be an ultra-high-finesse 
macroscopic resonator in the optical domain. Although the expected frequency shift for the 
lowest-order modes of such a resonator is near the resolution limit and is likely to be swamped 
by mirror imperfections, the quadratic dependence of the correction term on the mode indices J 
and V should allow the correction term to be relatively easily observed for higher-order modes. 
Verification of the present results in such a system would constitute a precision test of diffrac- 
tion and propagation in a resonator. 
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